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Abstract 


The  weight  function  procedure  allows  one  to  convert 
stress  intensity  factors  K  and  crack  displacement  informa¬ 
tion  obtained  for  one  crack  configuration  and  loading  into 
the  stress  intensity  factor  solution  for  the  same  geometry 
and  another  loading. 

The  feasibility  of  using  the  weight  function  idea 
for  a  two-dimensional  case  with  experimental  results  is 
demonstrated  in  this  work.  Mode  I  stress  intensity  factor 
K-j.  measurements  obtained  by  a  laser  interferometric  tech¬ 
nique,  and  "crack  mouth"  opening  displacement  measurements 
were  taken  for  an  edge  cracked  strip  subjected  to  four- 
point  bending.  These  results  were  used  to  construct  (numeri¬ 
cally)  a  weight  function  with  the  aid  of  a  computer  program 
written  for  this  purpose. 

Results  of  Kj  for  the  same  geometry  with  two  differ¬ 
ent  loading  configurations,  uniform  tension  and  three- point 
bending  (with  two  different  length  to  width  ratios)  were 
computed . 

These  results  agree  favorably  with  the  known  solu¬ 
tions  and  demonstrate  that  a  set  of  experiments  for  a  single 
loading  can  accurately  predict  the  stress  intensity  factor 
for  any  other  loading  configuration  of  the  same  geometry. 

The  advantage  of  the  weight  function  method  would  be  par¬ 
ticularly  important  if  these  loading  configurations  are 
difficult  or  impossible  to  reproduce  in  the  laboratory. 
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AN  EXPERIMENTAL  WEIGHT  FUNCTION  METHOD  FOR 
STRESS  INTENSITY  FACTOR  CALIBRATION 


I.  Introduction 


Background 

Fracture  mechanics  presently  provides  the  best 
available  tools  to  quantitatively  assess  the  influence  of 
preexistent  cracks  in  structures.  Although  fracture 
mechanics,  as  we  now  know  it,  is  a  relatively  new  tech¬ 
nology  (since  the  mid  1950s) ,  the  basic  ideas  were  already 
presented  by  Griffith  (Ref  1)  in  the  early  1920s. 

Current  linear  elastic  fracture  mechanics  (LEFM) 
concepts  assume  that  the  stress  intensity  factor  K,  the 
parameter  that  relates  load,  crack  length  and  geometry  con¬ 
trols  fracture  (K  =  =  constant  at  fracture)  and  crack 

propagation  (Ref  23) . 

Stress  intensity  factor  calibrations  are  required 
for  any  fracture  or  crack  propagation  analysis.  Since 
crack  tip  analysis  can  be  complex,  experimental  K  solutions 
often  become  necessary  to  verify  or  supplement  analytical 
or  numerical  solutions. 

For  complex  type  geometries  and  loading  configura¬ 
tions  which  commonly  exist  in  aircraft  structures,  experi¬ 
mentation  may  be  the  primary  method  to  obtain  a  reasonably 
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accurate  and  dependable  solution.  Conducting  experiments 
on  actual  parts  and  loadings  is  often  a  very  difficult  and 
costly  procedure. 

The  main  motivation  for  this  thesis  investigation 
was  to  find  a  technique  that  may,  in  some  cases,  greatly 
reduce  the  experimental  effort  required  in  order  to  obtain 
a  K  solution. 

An  analytical  approach  to  determine  the  stress 
intensity  factors  K  was  discussed  by  Bueckner  (Ref  2)  and 
Rice  (Ref  3) .  They  showed  that  once  the  displacement  field 
and  stress  intensity  factors  are  known  for  one  geometry  and 
loading  configuration  referred  to  subsequently  as  Case  1,  K 
may  be  obtained  for  any  other  loading  (Case  2)  applied  to 
the  same  crack  geometry  (Refs  4-8,  21).  The  method  depends 
mainly  on  the  reciprocal  theorem  and  strain  energy  expressions. 
A  load  and  material  independent  weight  function  can  be  con¬ 
structed  to  relate  the  stress  intensity  factor,  material 
properties,  and  crack  length  together  with  crack  surface 
opening  (for  more  details  see  Part  II  and  Appendix  A) . 

The  significant  advantage  of  this  procedures  lies  in 
the  fact  that  only  the  first  problem  (Case  1)  needs  to  be 
solved  directly.  Construction  of  the  weight  function  for 
this  case  allows  ready  determination  of  K  for  any  other 
loading  on  the  same  crack  geometry  (Cases  2,  3,  etc.). 

Stress  intensity  factors  may  then  be  obtained  for  a  variety 
of  other  specimen  loadings  some  of  which  may  be  highly  com¬ 
plex  compared  to  the  original.  In  addition,  the 
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computational  expenses  are  minimal  because  only  a  single 
relatively  simple  load  configuration  needs  to  be  calcu¬ 
lated. 

The  only  additional  information  required  for  any 
other  case  is  the  stress  distribution  at  the  crack  loca¬ 
tion  for  a  noncracked  body  under  the  load  configuration  in 
question.  This  stress  distribution  is  available,  in  many 
cases,  from  the  static  crack  free  stress  analysis  or  can 
be  obtained  with  much  less  effort  than  a  cracked  body 
analysis  requires. 

Purpose 

Although  weight  functions  have  been  employed  with 
analytical,  finite  elements  and  other  numerical  methods, 
the  author  is  unaware  of  attempts  to  combine  weight  func¬ 
tions  with  experimental  procedures. 

The  objective  of  this  work  was  to  show  that  the 
weight  function  idea  can  assist  the  experimentalist  by 
allowing  him  an  opportunity  to  obtain  the  stress  intensity 
factor  solution  for  a  complex  loading  configuration  from 
the  results  of  a  simple  experiment,  and  also  to  obtain  K 
solutions  for  a  variety  of  loadings  from  a  single  set  of 
experiments.  This  will  greatly  reduce  the  complexity  and 
scope  of  the  experimental  effort. 

General  Approach 

The  weight  function  procedure  is  demonstrated  with 
experimental  results  for  an  edge  crack  specimen. 
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The  reference  (Case  1)  problem  has  been  chosen  as 


an  edge  cracked  strip  subjected  to  four-point  bending. 

Crack  surface  displacement  measurements  near  the  crack  tip 
were  obtained  by  a  laser  interferometric  method  following 
the  technique  developed  at  the  Air  Force  Materials  Labora¬ 
tory  (Refs  9,  10,  19).  These  results  gave  the  mode  I 
stress  intensity  factor  K^.  for  that  case.  Crack  mouth 
opening  measurements,  using  a  clip  gage,  were  then  used 
to  construct  the  crack  surface  displacement  function 
incorporating  Orange  (Ref  11)  conic  section  approximation. 
These  experimental  results  were  used  to  evaluate  the  weight 
function  and  predict  the  stress  intensity  factor  Kj.  for 
other  cases,  which  were  chosen  to  be  uniform  tension  and 
three-point  bending.  The  results  were  compared  with  the 
known  solution  given  in  (Refs  12-14). 
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II.  Theo: 


General  Assumptions 

1.  Linear  elastic  fracture  mechanics  (LEFM)  applies 

i.e. ,  the  amount  of  placticity  near  the  crack  tip  is  rela¬ 
tively  small 


fe  <  o.i 

cl 

where  r^  is  the  radius  of  plastic  zone.  (See  Hertzberg, 

Ref  15,  and  Appendix  B  for  a  more  detailed  discussion.) 

2.  The  problem  is  assumed  to  be  purely  two- 
dimensional;  i.e.,  no  variations  through  the  thickness  are 
considered. 

3.  Body  forces  are  assumed  to  be  negligible. 

4.  Only  pure  Mode  I  crack  openings  will  be  dis¬ 
cussed. 

Stress  Intensity  Factor 

The  stress  intensity  factor  K  is  the  linear 
elastic  fracture  mechanic  (LEFM)  parameter  that  relates 
load,  crack  length  and  geometry.  The  LEFM  approach  to  pre¬ 
dicting  crack  growth  assumes  that  K  controls: 

1.  Fracture  (K  =  K^,  =  constant  at  fracture)  . 

2.  Crack  propagation  due  to  fatigue. 

3.  Crack  propagation  due  to  stress  corrosion. 
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Three  inodes  of  stress  intensity  factors  are  defined  depend¬ 
ing  on  the  crack  tip  opening  mode  (Fig.  2.1).  They  may  be 
stated  as: 


Mode 

I 

Opening 

KI 

=  lim 
r+0 

/JtFt 

a 

y 

(0=0) 

(2-1) 

Mode 

II 

Sliding 

H 

H 

K 

n 

*1  h- 
o  3 

/ Tnr 

axy 

(0=0) 

(2-2) 

Mode 

III 

Tearing 

KIII 

=  lim 
r-+0 

Siirr 

ayz 

(0=0) 

(2-3) 

where  crack  tip  coordinates  are  defined  in  Fig.  2.2 


MODE  I 
(Opening) 


MODE  II 
(Sliding) 


MODE  III 
(Tearing) 


Fig.  2.1.  Crack  Opening  Modes 
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Fig.  2.2.  Crack  Tip  Stress  Coordinates 

This  work  deals  only  with  the  opening  Mode  I, 
though  it  can  be  easily  extended  to  other  "modes." 


The  Weight  Function 

The  Mode  I  stress  intensity  factor  Kj  under  arbi¬ 
trary  loading,  shown  in  Fig.  2.3,  can  be  given  as: 


(2-4) 


(see  Appendix  A  and  Refs  2-6  for  detailed  discussion) . 
where : 

T  is  any  path  chosen  around  the  body  that  includes 
all  the  surface  tractions  and  body  forces, 

A  is  the  region  defined  by  T, 

T  designates  surface  tractions. 
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T  (s) 


Fig.  2.3.  Cracked  Body  Loading  Configuration 


Kj  is  the  Mode  I  stress  intensity  solution  corres¬ 
ponding  to  the  loads  T  and  f, 

f  are  body  forces. 


h  =  h(x,y,a)  =  (the  weight  function) 

1  (2-5) 


n(x,y,a)  is  the  c rack  vertical  displacement  for  any 
chosen  load  configuration  (independent  of  T 
and  f) , 


Kj.  in  the  weight  function,  Eq  (2-5)  ,  is  the  Mode  I 
stress  intensity  solution  that  corresponds  to 
the  n  solution ,  and 


for  plane  stress 
for  plane  strain. 


The  weight  function  h  is  independent  of  the  loads  T,  f  and 
material  properties  and  can  be  calculated  for  the  same 
geometry  from  another  known  solution  of  Kj  and  3n/3a. 

Using  the  principle  of  superposition  of  linear 
elasticity,  the  calculations  of  the  stress  intensity  factor 
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considering  the  case  of  a  cracked  body  under  a  certain 
loading  configuration  is  equal  to  the  case  of  the  cracked 
body  where  only  the  crack  surface  is  subjected  to  tractions 
p(x)  (coordinates  are  chosen  so  that  the  crack  is  in  x 
direction).  Here  p(x)  is  the  stress  that  will  occur 
at  the  crack  location  for  an  uncracked  body  subjected  to 
the  original  loading  configuration.  This  idea  is  demon¬ 
strated  in  Fig.  2.4. 


Fig.  2.4.  Superposition 

K=  K  +  K_ 

r®  r®  z© 


Part  B  can  be  viewed  as  a  noncracked  body;  therefore. 


If  one  chooses  the  path  r  as  any  path  that  includes 
the  crack  surface  and,  for  simplicity,  takes  a  straight 
crack  in  the  x-direction  as  in  Fig.  2.4  (the  last  assump¬ 
tion  is  not  essential  and  is  used  only  to  simplify  the 
equations) ,  then  the  Eq  (2-4)  for  K^.  can  be  stated  as: 

a 

h(x,a)p(x)dx  - 

Reversing  the  integration  limits  and  adding  the  two  inte¬ 
grals  we  get: 


f 


h(x,a)p(x)dx 


(2-6) 


K 


I 


h(x,a)p(x)dx 


(2-7) 


The  unknown  stress  intensity  solution  for  the  load¬ 
ing  p(x)  is  considered  Case  2,  and  the  known  solution,  from 
which  the  weight  function  will  be  computed,  is  taken  as 
Case  1,  then  Eq  (2-7)  can  be  written  as 


K 


I 


CD 


(2-8) 


and  from  Eq  (2-5) 


(2-9) 


If  Eq  (2-9)  is  substituted  into  Eq  (2-8) ,  the  result  is: 


(2-10) 


where 


K_  is  the  Mode  I  stress  intensity  factor  to  be 
found, 

P(x)q  is  the  stress  at  the  crack  location  for  a 

noncracked  body  with  loading  configuration  of 
Case  2, 


is  the  Mode  I  stress  intensity  factor  for 
Case  1  and  will  be  determined  experimentally, 
and 


is  the  derivative  with  respect  to  crack 
length  of  the  crack  opening  profile  for 
Case  1.  n  is  also  determined  experimentally. 


K  Calibration  by  Crack 
Tip  Displacement 

The  crack  tip  displacement  field  is  defined  by  the 
following  relations  that  represent  the  elasticity  solution 
for  the  stress-displacement  field  in  terms  of  K  near  the 
crack  tip  of  a  cracked  body.  One  can  take  these  displace¬ 
ment  expressions  for  Mode  I  opening  (Refs  17,  18,  10)  as: 


u 


u 


x 


y 


cos 


sin 


0  Tl-V  .  .2 

2  !_■ 1+v  +  sin 

0  r  2  2 

2  |l+v  "  COS 


(2-11) 


(2-12) 


for  plane  stress;  and 
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\  *  TT  VT  I  [1-2v  +  si"‘  f]  <2-13> 

uy=_^Vsr  sin  !  [2-2v  -  cos!  |J  (2-14) 


for  plane  strain. 

Consider  the  vertical  displacement  u^  along  the 
crack  surface  for  which 


Q  =  TT 
r  =  t 


where  t  is  a  coordinate  along  the  crack,  starting  at  the 
crack  tip.  One  should  note  to  conform  with  usual  practice, 
n  replaces  the  displacement  in  y  direction  u^,.  Recall  from 
elasticity  that: 


G  2 ( 1+v) 


(2-15) 


Substituting  into  Eqs  (2-12)  and  (2-14)  the  crack  displace¬ 
ment  for  plane  stress  and  plane  strain  can  be  written  as 


(2-16) 

(2-17) 


The  equation  holds  only  very  near  the  crack  tip 


t  <<  a 
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and  allows  one  to  calibrate  by  measuring  crack  opening 
at  a  known  distance  very  near  the  crack  tip. 

Crack  Opening  Calibration 

In  order  to  find  the  slope  with  respect  to  the 
crack  length  3n/3a  of  the  crack  opening  profile,  an  approxi 
mation  has  been  used.  Orange  (Ref  11)  suggested  a  conic 
section  approximation  for  a  finite  width  edge  cracked 
plate.  This  approximation  is  based  on  assuming  a  general 
conic  section  function  for  the  crack  opening  shape  and 
fitting  the  crack  tip  radius  of  curvature  and  crack  mouth 
opening  to  determine  the  unknown  coefficients.  Further 
development  of  this  relationship  is  presented  in  Appendix  E 
Following  Orange  and  substituting  for  the  crack  tip  curva¬ 
ture  in  terms  of  the  stress  intensity  factor,  using  crack 
tip  displacement  equations,  the  following  expression  is 
derived: 


(2-18) 


where  nQ  is  the  crack  mouth  displacement  and  m  is  the  conic 
section  coefficient. 


m  =  it 


-2 


(2-19) 


Y  is  a  nondimensional  stress  intensity  factor 


a/ a- 


(2-20) 
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Fig.  2-5.  Edge  Crack 


Crack  displacements  of  finite  single  edged  cracked 
strips  under  bending  or  tension  generated  by  the  Orange 
Eq  (2-18)  were  compared  to  the  collocation  method  pre¬ 
sented  by  Gross  (Ref  20)  and  found  to  match  within  3.2  per¬ 
cent.  The  Orange  equation  was  compared  to  the  Rice  (Ref  3) 
solution  for  a  center  crack  considering  an  infinite  plate. 
The  results  were  found  to  be  exactly  the  same.  The  Orange 
equation  was  also  used  successfully  by  Grandt  (Refs  6,  21) 
for  cracks  in  fastener  holes  and  for  radially  cracked  rings. 
Substituting  the  expression  for  Y  Eq  (2-20) ,  and  m  Eq 
(2-19)  into  Eq  (2-18) ,  we  obtain 


(2-21) 


(2-22) 
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In  order  to  simplify  the  above  equation. 


let 


and 


a2  7TH2a 


Thus, 


n  =  [Bt  +  Rt2]' 


Also, 


In 

8a 


1 

2 


-h 

[Bt  +  RT2] 


h  1Bt  +  Rt!> 


(2-23) 


(2-24) 


(2-25) 


(2-26) 


From  Fig.  2.5,  if  t 


a-x. 


then  9a  1* 


Thus, 


in  =  I  rBt 

9a  2  [Bt 


+  Rt2] 


-Js 


.  9B 
9a 


t  +  B  + 


|5  t’  + 

9a 


2Rt]  (2-27) 


where: 


16Kj 

ttH2 


(2-28) 


and 


9R 

"V  „ 

9^  2n  2 

o  o 

16K^ 

9Kj. 

9a 

ttH  2  a2 

a2 

9a  3 

a 

irH2a 

9a 

(2-29) 


Eq  (2-27)  therefore  defines  9n/9a  in  terms  of  the  crack 
mouth  displacement  no  and  the  stress  intensity  factor  Kj 
and  their  derivative  with  respect  to  the  crack  length. 


Kj  For  Case  2 

Substituting  Eq  (2-27)  into  Eq  (2-10)  and  also 
changing  the  integration  variable  from  x  to  t  one  obtains 


K, 


H 


[© 


c© 


*  r\ 


[Bt+Rt2]  +(|f  +  2R)t +  |^  t2]dt  (2-30) 

da  d  cl 


9B  9  R 

where  B,  R,  r-  and  -5 —  are  defined  previously  and  calcu- 

da  0  3 


lated  for  Case  1. 

From  Eq  (2-30)  it  is  clear  that  measuring  the 
stress  intensity  factor  Kj  and  the  crack  mouth  opening  riQ 
for  various  crack  lengths  and  calculating  their  derivatives 
with  respect  to "a"  allows  one  to  construct  the  weight  func¬ 
tion  for  the  Case  1  loading  configuration.  By  knowing 
the  stress  distribution  at  the  crack  location  for  a  non- 
cracked  body  under  the  loading  configuration  of  Case  2 
p,tb  ,  one  may  obtain  for  Case  2. 

Computing  this  integral  is  quite  tedious  though 
possible  for  very  simple  forms  of  p(t) .  However,  numerical 
integration  is  possible  for  any  p(t)  even  if  represented 
by  a  discrete  numerical  solution.  In  performing  numerical 
integration,  one  should  be  cautioned  to  the  fact  that  the 
weight  function  and  therefore  the  whole  expression  under  the 
integral  is  square  root  singular  at  the  crack  tip  t=0  but 
still  it  can  be  shown  that  for  the  problems  studied  here 
a  limit  exists  for  the  integral  and,  in  most  cases,  the 
numerical  integration  converges  quite  rapidly. 
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III.  Experimental  Technique 

The  Laser  Interferometry  Method 

Stress  intensity  factor  K  calibrations  obtained 
by  measuring  crack  tip  displacements  require  a  technique 
capable  of  accurately  measuring  extremely  small  displace¬ 
ments  very  near  the  crack  tip.  As  a  "rule  of  thumb" 
the  relations  for  the  crack  tip  displacement  field  as 
expressed  in  Eqs  (2-11  to  2-14)  are  regarded  to  be  suf¬ 
ficiently  accurate  within  a  distance  of  a/20  or  smaller 
from  the  crack  tip,  where 'a*' is  the  crack  length.  The  laser 
interferometry  technique  described  here  was  reported  in 
Refs  9,  10,  and  19  to  measure  displacements  of  0.1  microns 
at  a  distance  of  50  microns  from  the  crack  tip. 

Results  were  reported  leading  to  values  that 
varied  up  to  15  percent  from  the  theoretical  value  (for 
individual  points) .  Most  of  these  deviations  were  attri¬ 
buted  to  the  imperfect  crack  and  material  characteristics 
and  not  to  measurement  errors.  The  technique  is  basically 
similar  to  conventional  interferometry  but  only  measures 
in-plane  displacements  on  the  surface. 

Two  small  indentations  are  placed  with  a  diamond 
indentor  on  both  sides  of  the  crack  near  the  crack  tip. 
Those  indentations  are  typically  square  based  pyramids 
with  a  base  length  of  20-40  microns  and  located  50  microns 
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apart  on  both  sides  of  the  crack  (see  Fig.  3.1). 


Fig.  3.1.  Surface  Indentation  (x280  magnification) 

A  laser  coherent  source  impinges  upon  the  indentations. 

The  beam  is  diffracted  back  at  an  angle  a q  with  respect  to 
the  incident  beam  as  shown  schematically  in  Fig.  3.2. 

Since  the  indentations  are  placed  close  together, 
the  respective  diffraction  beam  overlaps  resulting  in  inter¬ 
ference  fringe  patterns  on  either  side  of  the  incident 
laser  beam.  A  photograph  of  a  typical  interference  pat¬ 
tern  is  shown  in  Fig.  3.3. 

The  relation  between  the  indentation  spacing  and 
the  fringe  order  is  set  by  the  following  equation  (Ref  9)  : 

d  sin  a  =  mX  (3-1) 
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where  d  is  the  distance  between  indentations,  X  is  the  wave 
length  of  the  source  and  m  is  the  fringe  order. 

As  a  load  is  applied  to  the  specimen,  the  crack 
opens  and  the  distance  between  the  indentation  changes. 

The  relation  between  the  change  in  d  and  the  change  in  the 
fringe  order  <$m  from  a  fixed  point  of  view  is 


fid  =  fim  X/sin  a 

o 


(3-2) 


Now  observing  the  fringe  pattern  from  a  fixed  point  allows 
one  to  count  the  number  of  fringes  passing  and  hence  deter¬ 
mine  the  change  in  the  fringe  order.  Eq  (3-2)  will  then 
give  the  change  in  the  distance  d  between  the  indentations 
that  correspond  to  the  crack  opening.  Averaging  the  left 
and  right  fringe  number  will  eliminate  free  body  motion. 


Test  Apparatus 

The  test  specimens  used  were  flat  strips  of 
Aluminum  7075  T651.  A  small  V  notch  was  placed  on  the  edge 
of  the  specimen  to  allow  precracking  and  a  shallow  groove 
was  cut  at  the  edge  as  to  fit  the  clip  gage  used  to  measure 
crack  mouth  openings.  (See  Fig.  3.4.) 


Fig.  3.4.  Test  Specimen 
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The  specimen  was  placed  in  a  four-point  bending  fixture. 
Load  was  applied  by  a  standard  Instron  machine  equipped 
with  a  compression  load  cell  (Fig.  3.5). 


Fig.  3.5.  Specimen  Set  Up  in  Bending  Fixture 

The  crack  mouth  opening  was  measured  by  a  "clip 
gage"  placed  in  the  premachined  groove.  The  "clip  gage" 
consists  of  two  small  steel  plates  on  which  strain  gages 
were  bonded  (Fig.  3.6).  Precalibrating  the  "clip  gage" 
allows  one  to  relate  the  mouth  opening  with  the  changes  in 
strain  measured. 

Two  pyramid  type  indentations  have  been  placed 
near  the  crack  tip  by  a  standard  Lietz  microhardener  tester 
with  a  diamond  indentor  of  a  square  based  pyramid  shape 
with  face  angles  of  136  degrees.  The  indentations  were 
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Fig.  3.6.  Clip  Gage  Setup 

impinged  by  a  Spectra  Physics  Model  120  5  MW  HeNe  laser. 

The  laser  beam  wave  length  X  was  0.6328  microns  and  had  a 
divergence  angle  of  0.71  mr  and  a  beam  diameter  of  0.57  mm. 
The  laser  was  located  so  that  the  beam  impinged  perpendicu¬ 
lar  to  the  specimen  surface.  The  angle  of  the  diffracted 
beam  aQ  can  be  calculated  directly  by  knowing  the  face 
angle  of  the  indentation.  (See  Fig.  3.7.) 

The  fringes  created  by  the  interference  of  the 
reflected  beams  were  picked  up  by  two  photoresistors  that 
were  located  in  a  fixed  place  on  the  fringe  pattern. 

The  active  face  of  the  photo  resistors  was  masked 
leaving  a  narrow  slit  (smaller  than  the  fringe  spacing) 
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As  load  was  applied  and  the  fringes  started  moving  the 
photoresistor  created  an  electrical  signal  that  after 
proper  amplification  could  be  recorded  with  the  aid  of  a 
strip  chart  recorder. 

Test  Procedure 

Two  Aluminum  7075  T651  specimens  were  each  used  to 
test  5-10  different  crack  lengths.  The  first  specimen  was 
used  for  longer  cracks  and  the  second  for  shorter  cracks 
with  some  overlapping  between  the  regions  that  allowed  com¬ 
paring  between  the  two.  The  specimens  were  precracked  in 
the  long  transverse  (LT)  direction  of  the  material  in  three- 
point  bending  on  a  Schenck  fatigue  machine  using  cyclic 
loads  at  a  frequency  of  20  Hz.  The  precracking  stress  inten¬ 
sity  factor  varied  from  6000-9000  psi  /In  .  The  higher  K^. 
was  used  to  initiate  and  propagate  the  crack  to  about  80 
percent  of  the  desired  length,  and  the  lower  K^.  levels  were 
applied  to  get  the  final  sharp  crack. 

Two  indentations  were  placed  on  the  surface  50-200 
microns  behind  the  crack  tip  as  explained  previously.  A 
load  of  200  grams  on  the  indentor  diamond  head  was  found 
most  suitable  to  create  good  reflecting  indentations  with 
a  base  square  dimension  of  approximately  30  microns.  The 
specimen  was  mounted  on  the  loading  fixture,  the  two 
photoresistor  output  signals  together  with  the  clip  gage 
and  load  cell  signals  were  monitored  by  a  four-channel 
strip  chart  recorder  to  create  a  common  time  basis. 
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The  specimens  were  loaded  typically  to  values  of 
Kj.  =  12-14  ksi  /In  so  the  amount  of  plasticity  in  the 
specimen  was  kept  small  (see  Appendix  B) .  Two  load  cycles 
were  applied  prior  to  actual  testing  so  that  the  crack 
would  appear  to  its  full  length  bn  the  surface  and  would 
easily  be  seen  under  a  microscope. 

During  the  preload  cycle  the  photoresistors  were 
relocated  to  a  position  such  that  the  maximum  fringe  inten¬ 
sity  could  be  observed.  Since  the  photoresistors  were 
fixed  in  position  while  the  specimen  with  the  indentation 
was  displacing  due  to  the  applied  load,  it  became  neces¬ 
sary  to  establish  the  final  location  of  the  resistors  under 
a  mean  load  so  that  the  indentations  could  still  stay 
within  the  beam  boundaries  over  the  extreme  deflection. 

Three  load  cycles  were  applied  to  each  specimen  for 
each  crack  length.  After  each  test  the  crack  length  was 
measured  in  order  to  account  for  any  crack  propagation  dur¬ 
ing  the  test  cycle.  Then  the  -rack  was  propagated  to  the 
next  stage  in  the  Schenck  fatigue  machine.  In  addition, 
for  two  different  crack  lengths  the  crack  surface  displace¬ 
ment  was  measured  along  the  crack  at  2-6  different  locations 
using  previously  applied  indentations  in  order  to  verify 
the  Orange  (Ref  11)  conic  section  approximation  for  the 
crack  opening. 
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IV.  Data  Reduction  and  Numerical  Technique 

Kj  For  Case  1^ 

The  fringe  patterns,  crack  mouth  opening  and 
loads  were  all  recorded  on  a  four-channel  strip  chart 
recorder.  A  typical  chart  for  a  single  loading  cycle  is 
shown  in  Appendix  C,  Fig.  C.l. 

As  explained  previously,  two  sets  of  fringes  were 
created  by  the  reflected  beam  from  the  inaentation.  For 
each  set  of  fringes,  a  slope  of  number  of  fringes  per  unit 
load  6m/p  can  be  obtained.  These  slopes  were  averaged  to 
eliminate  free  body  motion.  A  more  detailed  description 
of  the  laser  interferometer  data  reduction  procedure  can 
be  found  in  Appendix  C. 

Using  Eq  (3-2) 

n  _  1  6d  _  1.  _  3.  6m  .  1  ,  .  . . 

P  2  P  P  2  p  sin  a  1 

o 

where  X/sin  a  is  the  calibration  factor  of  the  inter- 
o 

ferometer  used  here.  Eq  (4-1)  allows  one  to  compute  the 
displacement  per  unit  load.  The  loading  and  unloading 
results  of  three  loading  cycles  were  averaged  to  give  an 
average  displacement  per  unit  load.  The  distance  t  of  the 
indentation  behind  the  crack  tip  was  computed  by  averaging 
the  initial  distance  before  the  test  with  the  one  measured 
after  three  test  cycles.  Since  only  three  cycles  were 
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applied  and  the  amount  of  crack  propagation  was  very  small 
(a  few  microns)  linear  interpolation  was  adequate  to  give  a 
sufficiently  accurate  estimate. 

Kj  for  Case  1  can  now  be  determined  from  the  crack 
displacement  near  the  crack  tip  Eqs  (2-12)  or  (2-14) .  In 
most  cases  0  is  close  enough  to  180°  so  the  simplified  form 
of  Eq  (2-17)  can  be  used. 

Since  the  measurements  are  on  the  surface,  the 
stress  field  should  be  considered  as  plane  stress  and  H 
should  be  given  the  value  of  E  (1x10^  for  aluminum  7075T651) . 
Now  Eq  (2-17)  takes  the  form 

tt  =  (?) !  (¥)" 

This  data  will  usually  be  converted  to  some  nondimensional 
form. 

Crack  Mouth  Opening 

The  crack  mouth  opening  h0  was  recorded  directly 
from  the  "clip  gage"  readings.  The  crack  mouth  opening  per 
unit  load  was  then  readily  obtained  from  a  nQ  versus  P  plot 
using  only  the  linear  part  of  the  curve  and  neglecting  the 
initial  nonlinear  part  which  is  associated  with  crack 
closure  effects. 

Numerical  Techn ique  for  K 
Case  2  Computation 

The  rest  of  the  computation,  which  involved  calcu¬ 
lating  stress  intensity  factors  for  different  Case  2 


27 


loadings  were  programmed  and  handled  numerically.  The 

Kt  data  was  fitted  by  a  suitable  interpolation  curve. 

© 

Interpolated  values  of  KT 

(D 

any  desired  crack  length  a/w.  In  a  similar  way  the  crack 
mouth  opening  data  for  Case  1  n  was  fitted  by  an  inter- 

polating  curve  that  allows  one  to  compute  values  of  n  and 
9no 

-r—  for  the  desired  points.  These  interpolations  are  calcu- 

0  a 

lated  in  a  different  subroutine  so  any  desired  interpola¬ 
tion  function  can  be  tried  without  changing  the  main  pro¬ 
gram.  A  comparison  table  for  Kj  and  nD  for  interpolated 
data  versus  actual  data  is  printed  out  so  one  can  judge  the 
accuracy  of  the  interpolation  process  chosen.  Fourth  order 
least  square  polynomials  will  give  reasonable  results  in 
many  cases. 

Using  Orange  (Ref  11)  conic  equation  discussed  in 

Appendix  E  for  the  crack  opening  profile  the  derivative  of 

the  crack  opening  along  the  crack  with  respect  to  the  crack 

length  can  be  computed  through  Eq  (2-27) . 

Since  different  factors  are  used  to  nondimen- 

tionalize  the  K  and  n  data,  care  should  be  taken  to  bring 
KI  °  9no 

all  the  Kj,  -^a  ,  and  -^a  ■  to  a  common  basis.  Using 

Eq  (2-9) ,  the  weight  function  is  computed  for  each  crack 
length  at  any  point  along  the  crack.  The  constant  H  is 
determined  depending  on  plane  stress  or  plane  strain  con¬ 
ditions  (see  Appendix  B) .  The  stress  distribution  along 
the  crack  location  for  a  noncracked  body  p(t)  is  supplied 


3Kt 

_i© 


were  computed  for 
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by  a  separate  subroutine  and  can  be  easily  fitted  with 
available  data  or  some  known  analytical  solution.  Kj.  for 
Case  2  now  is  computed  through  integrating  Eq  (2-30). 

The  integral  is  performed  numerically  using  a  Rom¬ 
berg  integration  scheme  (Ref  24) I  Although  the  weight 
function  and  thus  the  integrand  in  Eq  (2-30)  is  singular  at 
the  crack  tip,  the  integral  will  still  converge  to  a  dis¬ 
tinct  value  quite  rapidly  for  most  p(t)  functions.  Because 
of  the  crack  tip  singularity  the  lower  limit  of  this  inte¬ 
gral  cannot  be  set  to  zero.  A  value  close  to  zero  should 
be  taken  and  this  value  can  be  made  smaller  and  smaller  by 
an  iterative  process  until  the  value  of  the  integral  con¬ 
verges  within  a  specified  accuracy. 

Since  the  numerical  integration  process  may  be 
executed  many  times  until  convergence  is  obtained,  an  effi¬ 
cient  algorithm  should  be  used  to  reduce  the  required  com¬ 
puter  time.  Romberg  integration  scheme  used  here  is  a 
powerful  and  efficient  numerical  integration  technique. 

It  is  based  upon  the  trapezoidal  rule  combined  with  Richard¬ 
son  extrapolation.  More  details  of  this  scheme  can  be 
found  in  Hornbeck  (Ref  24) .  Using  this  scheme  reduced  the 
computer  time  by  a  factor  of  about  200  compared  to  the 
Simpson's  rule  integration  scheme.  A  detailed  description 
of  the  computer  program  is  presented  in  Appendix  D. 
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V.  Experimental  Results 

Experimental  Measurements 

Two  Aluminum  7075  T651  edge  cracked  strips  were 
used  in  the  experiments.  The  specimens  were  cut  out  from 
a  rolled  plate  and  the  era  ’  was  grown  in  the  long  transverse 
(LT) direction  of  the  material. 

The  dimensions  of  the  specimens  and  the  region  of 
crack  length  in  which  they  were  tested  are  shown  in  Fig. 

5.1. 
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N* 

Q 

F 

a 

Specimen  A 

9" 

0.24" 

1" 

.27” 

0.2328"  -0.7517" 

Specimen  B 

9" 

0.245” 

0.7885" 

.31” 

0.036” 

0.036” 

0 . 08153"-0 .3753" 

*N  dimension  determined  for  clip  gage  purposes. 


Fig.  5.1.  Specimen  Dimensions 
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Specimen  A  was  tested  with  ten  different  crack 
lengths  with  crack  length  to  width  ratios  ranging  from 
.2328-. 7517.  Specimen  B  was  tested  with  five  different 
crack  lengths  where  a/w  ranged  from  .1034-. 476.  Both 
specimens  were  loaded  in  a  four-point  bending  fixture  as 
shown  in  Fig.  5.2. 


Fig.  5.2.  Loading  Fixture  Dimension 


Each  specimen  was  loaded  through  three  loading  cycles,  and 
near  crack  tip  displacements  (by  laser  interferometry) 
and  crack  mouth  openings  were  recorded.  Table  I  represents 
the  average  results  for  all  crack  lengths  measured. 


Data  Interpolation 

The  stress  intensity  data  calculated  from  the 
interferometry  measurements  was  converted  to  a  nondimen- 
sional  form  of 


Kj  KjbW2 

o/ira  6M/5fa 


(5-1) 
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TABLE  I 


Experimental  Measurement  Results 


Specimen 

Nondimensional 
crack  length 

Distance  of 
indentation 
from  crack 
tip 

No  of  fringes 
per  unit  load 

Crack  mouth 
displacement 
per  unit  load 

a/w 

r 

6m/p 

n0/P 

Microns 

ms 

B 

.1034 

101 

5.093 

0.568 

B 

.1586 

83.5 

6.623 

0.7453 

B 

.2064 

88 

7.6713 

0.9836 

A 

.2328 

121.5 

7.3128 

0.856 

B 

.2646 

86.5 

9.116 

1.268 

A 

.2666 

99.5 

7.580 

1.14 

A 

.3094 

143 

10.180 

1.405 

A 

.346 

56 

7.340 

1.705 

A 

.3978 

72 

9.140 

2.358 

A 

.436 

46.5 

8.0146 

2.727 

B 

.476 

47 

12.808 

4.121 

A 

.4893 

50 

10.130 

3.6585 

A 

.5780 

111.5 

18.354 

6.04 

A 

.6467 

74 

20.460 

9.216 

A 

.7517 

70 

35.500 

20.681 

p  •  1 

where  M  is  the  bending  moment  M  =  — ^ —  . 

The  crack  mouth  opening  data  was  converted  to  a  nondimen- 
sional  form 


o 

2aa 


H  = 


n  bW  H 
o _ 

12Ma 


(5-2) 


The  nondimensional  data  for  K^.  and  no  was  interpolated  by 
a  least  square  fourth  order  polynomial.  Table  II  and 
Table  III  present  the  actual  experimental  and  the  interpo¬ 
lated  data  for  Kj  (nondimensional)  and  no  (nondimensional) 
respectively  including  the  relative  difference  in  percent 
of  the  interpolating  polynomial  at  the  test  point. 

Fig.  5.3  presents  the  actual  experimental  nondimen¬ 
sional  Kj  compared  to  the  known  reference  results  from 
Tada  (Ref  12).  Fig.  5.4  presents  the  actual  experimental 
nondimensional  data  for  riQ  compared  to  the  known  reference 
results  from  (Ref  12) .  The  average  difference  between  the 
reference  values  of  Kj  and  the  experimental  data  is  6%  while 
the  maximum  difference  for  a  single  measurement  is  up  to 
14%.  These  results  are  similar  in  accuracy  to  the  ones 
obtained  by  Macha  (Ref  10)  for  the  laser  interferometry 
technique.  The  values  of  Kj  from  the  interpolating  poly¬ 
nomial  that  were  actually  used  to  compute  the  weight  func¬ 
tion  are  similar  in  accuracy  with  a  maximum  difference  of 
13%  and  an  average  difference  of  6%  compared  to  Ref  12. 

One  should  note  that  the  relatively  large  difference  in 
the  slope  of  -r— -  for  some  of  the  short  crack  lengths 


TABLE  II 


1 


K_  (Nondimens iona I)  Comparison  of  Actual  Data  to 
Interpolated  Data 


Actual 

Interpolated 

Relative 

Difference 

a/W 

KjbW  2 

KTbW2 

* 

in  % 

oM/ira 

aM/ira 

.1034 

.  9073 

.  9467 

4.34 

.1586 

1.0575 

.  9776 

-7.56 

.2064 

1.0572 

1.0545 

-.26 

.  2328 

1.1137 

1.1060 

-.69 

.2646 

1.1007 

1.1705 

6.35 

.2666 

1.1862 

1.1746 

-.98 

.3096 

1.2366 

1.2590 

1.81 

.3460 

1.3490 

1.3244 

-1.82 

.3978 

1.3804 

1.4086 

2.04 

.4360 

1.4388 

1.4696 

2.14 

.4760 

1.5646 

1.5424 

-1.42 

.4893 

1.6551 

1.5710 

-5.08 

.5780 

1.8481 

1.8774 

1.59 

.6467 

2,3460 

2.3688 

.97 

.7517 

3.9513 

3.9413 

-.25 

♦The  interpolating  polynomial  is: 


— -  =  70.75(a/W)  4  -92.45(a/W)  3 +42. 95 (a/W) 2-6. 53 (a/W) +1.257 

aM/ira 


! 

! 


TABLE  III 


n  (Nondimensional)  Comparison  of  Actual  Data  to 
°  Interpolated  Data 


gg 

Actual  T) 

o 

Interpolated  nQ 

Relative 
Difference 
in  % 

a/W 

nQbW2H 

"obW’H  . 

12Ma 

12Ma 

.1034 

1.9440 

1.9500 

.31 

.1586 

1.6643 

1.6491 

-.91 

.2064 

1.6855 

1.6334 

-3.09 

.2328 

1.6163 

1.6879 

4.43 

.2646 

1.6950 

1.7933 

5.80 

.2666 

1.8796 

1.8011 

-4.18 

.3096 

1.9961 

1.9923 

-.19 

.3460 

2.1660 

2.1823 

.75 

.3978 

2.6055 

2.4945 

-4.26 

.4360 

2.7493 

.  2.7696 

.74 

.4760 

3.0632 

3.1251 

2.02 

.4893 

3 . 2866 

3.2648 

-.66 

.5780 

4.5933 

4.6431 

1.08 

.6467 

6.6244 

6.5725 

-.78 

.7517 

12.0980 

12.1082 

.08 

*The  interpolating  polynomial  is: 
n  bWH 

=  210.38 (a/W) 4 -2 8 2. 57 (a/W) 3 +150. 74 (a/W) 2 -32 . 16 (a/W) + 
+3.952 
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Kj  bWz 

6M/ia 

-  Tada  (Ref  12) 

-  Interpolation  curve 

for  the  experimental 
data 

A  Experimental  data. 
Specimen  A 

A  Experimental  data,  < 
Specimen  B  ' 

Fig.  5.3.  Stress  Intensity  Factor  (Nondimensional) 

for  Pure  Bending 
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(see  Fig.  5.3)  will  have  only  a  minor  influence  on  the 
weight  function  computation  as  can  be  seen  from  examining 
Eq  (2-30) .  nQ  data  as  well  as  the  interpolating  polynomial 
fits  the  reference  (12)  values  much  better  as  seen  in 
Fig.  5.4. 

Crack  Profile  Measurements 

For  two  different  crack  lengths  on  specimen  A  the 
crack  opening  was  measured  along  the  crack  at  several 
points  in  addition  to  the  point  near  the  crack  tip  and  at 
the  crack  mouth.  These  measurements  were  done  by  the  laser 
interferometry  technique.  The  existing  indentations  left 
behind  as  the  crack  propagated  were  used  for  these  measure¬ 
ments.  The  crack  openings  along  the  crack  were  compared  to 
the  Orange  (Ref  11)  conic  sections  and  found  to  match  this 
approximation  excellently  (see  Figures  5.5  and  5.6). 

Case  2  Loadings 

Three  different  loading  configurations  were  chosen 
for  Case  2,  uniform  tension,  three -point  bending  with  S/W 
of  4  and  8  (see  Fig.  5.7).  Each  of  the  different  Case  2 
loading  configurations  requires  another  function  p(t)  as 
the  noncracked  body  stress  at  the  crack  location  stresses. 

Uniform  Tension 

In  this  case, 

p(t)  =  const.  =  a 

o 

where  a  is  the  remote  stress.  Since  the  standard 
o 
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Fig.  5.5.  Nondimensional  displacement  along  Crack  Surface  Compared  to  Orange  Equation 

(Ref  11)  Pure  Bending  a/W=0.436 


Uniform  Tension 


i 


Three-point  Bending 


Fig.  5.7.  Case  2  Loading  Configurations 


nondimen sional  representation  for  for  uniform  tension  is 


K 

Tiaa 

o 

the  function  p(t)  should  be  factored  by  l//rr aaQ 
Kj  in  the  desired  form.  Therefore,  take 


(5-4) 

to  get 


P(t)  =  -i-  (5-5) 

✓SiT 

(in  the  computer  program  l/Zifa  is  already  factored  in  the 
weight  function  so  p(t)  is  substituted  as  1.). 

The  results  of  the  nondimensional  stress  intensity 
factor  for  uniform  tension  at  equally  spaced  points  from 
a/W=.l  to  a/W=.7  are  presented  in  Table  IV.  These  results 
are  compared  to  Tada's  (Ref  12)  solution.  The  results 
agree  with  the  reference  solution  within  8  percent,'. 
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TABLE  IV 


Kj.  (Nondimensional)  for  Uniform  Tension 


Crack  Length 

a/W 

KI 

KI 

talCUlaLcU 

a/rra 

Using  the 
Weight  Function 

a/ ira 

Ref  12 

Relative 
Difference 
in  % 

.1000 

1.2705 

1.1957 

6.26 

.1500 

1.2697 

1. 2682 

.12 

.2000 

1.4041 

1.3667 

2.74 

.2500 

1.5540 

1.4941 

4.01 

.3000 

1.6855 

1.6551 

1.83 

.3500 

1.8170 

1.8565 

-2. 13 

.4000 

1.9967 

2. 1080 

-5.28 

.4500 

2.2917 

2.4241 

-5.46 

.5000 

2.7719 

2.8266 

-1.  93 

.5500 

3.4752 

3.3486 

3.78 

.6000 

4.3636 

4.0432 

7.92 

.6500 

5.3158 

4.9993 

6.33 

7000 


5.1832 


6.3755 


-3.02 


Fig.  5.8  compares  the  reference  solution  (Ref  12) 


to  the  values  calculated  at  the  same  a/W  points  where 
experiments  were  performed  on  the  four-point  bending  speci¬ 
men.  One  can  see  that  the  last  point  is  the  only  calculated 
value  that  does  not  agree  favorably  with  the  reference  data. 
Since  the  computations  required  differentiation  of  experi¬ 
mental  data,  one  should  be  very  suspicious  about  the 
results  obtained  at  the  experimental  interval  end  points. 
Close  to  the  end  points  computing  derivatives  actually 
requires  extrapolations  that  can  be  very  inaccurate,  since 
the  behavior  of  the  function  is  unknown  in  these  regions. 

Investigating  Eq  (2-30)  shows  clearly  that  the 
derivatives  3K/3a  and  3nQ/3a  have  a  relatively  larger  influ¬ 
ence  on  the  value  of  the  weight  function  at  longer  crack 
lengths;  therefore,  the  computed  values  are  more  likely  to 
get  distorted  at  the  upper  limit  of  the  experimental  a/W 
interval  as  can  indeed  be  seen  in  Fig.  5.8. 

Three-Point  Bending 

A  solution  for  the  stress  distribution  of  a  non- 
cracked  strip  under  three-point  bending  in  a  cross-section 
through  the  center  is  given  by  Timoshenko  (Ref  27)  : 


(5-6) 
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The  terms  of  this  equation  are  defined  in  the  fol¬ 
lowing  Figure  5.9.  a  is  the  stress  in  x- direction  along 

X 

AD. 


p 

A 

y 

X 

2c=W 

D 

P/2 

2Jl=S 

Fig.  5.9.  A  Uniform  Strip  Under  Three  Point  Bending 

According  to  Timoshenko,  this  expression  gives  the  stress 
with  very  good  accuracy  except  for  point  D  where  an  error 

p 

of  .121  — r-  exists  based  on  a  more  accurate  solution.  This 
cb 

term  will  lead  in  our  worst  case  to  an  error  of  4  percent 
in  the  load  distribution  at  the  crack  mouth  location,  but 
since,  at  this  point  the  weight  function  has  its  lowest 
value,  the  integral  of  the  product  of  h  and  p(t)  is  expected 
to  have  much  better  accuracy.  The  standard  nondimensional 
solution  is  of  the  form 


K 


I 


a  /Ira 


(5-6a) 
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where  a  =  — —  -  and  M  =  —■ 
4  be  2  2 


So  the  undimensional  form  of  becomes 


Kj  4c2b 


/ira  3PJ l 


(5-7) 


The  factor  l/Zira  is  already  included  in  our  computer  pro¬ 
gram  in  the  weight  function  so  the  form  of  p(t)  that  is 
required  in  order  to  get  Kj  in  the  form  of  Eq  (5-7)  is 

a  (t)4bc2 

p(t)  =  -X  TpI -  (5_8) 


Substituting  Eq  (5-6)  to  Eq  (5-8)  yields 

I  *  l  (f)(i)'  l5'9» 

where  S  =  2£ 

and  W  =  2c 


p(y) 


=  (c 
3tt  \  £ 


Changing  to  crack  tip  coordinates  t  the  expression 
for  y/c  becomes 


y/c  = 


(W/2) 


-a)  =  2  (s  -  §) 


+  1 


(5-10) 


Substitute  Eq  (5-10)  into  Eq  (5-9)  and  get 


(5-11) 
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Results  for  the  nondimensional  stress  intensity 
factor  considering  three-point  bending  (S/W=4  and  S/W=8) 
are  presented  in  Tables  V  and  VI  respectively.  These 
results  are  compared  with  known  results  from  Ref  12  and 
Ref  28.  The  reference  solution  claims  to  be  accurate 
for  a/W  <  6  and  in  this  interval  the  calculated  experi¬ 
mental  values  are  within  6.15  percent  accuracy  for  S/W=4 
and  5.9  percent  for  S/W=8. 

Fig.  5.10  and  Fig.  5.11  show  the  calculated  K^. 
for  three-point  bending  for  the  a/W  points  where  actual 
data  was  gathered  for  the  four-point  bending.  As  for  the 
uniform  tension  case  the  last  point  is  the  only  one  not 
fitting  the  expected  curve  because  of  the  inaccurate 
extrapolation  procedure  explained  previously. 
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TABLE  V 


Kj  (Nondimensional)  for  Three-Point  Bending,  S/W=4 


Crack  Length 

a/W 

Kj  (Nondimensional) 

Calculated  Using 
the  Weight  Function 

Kj  (Nondimensional) 
Ref  28 

Relative 
Difference 
in  % 

.1000 

1.0134 

.9849 

2.89 

.1500 

.9447 

.9731 

-2.91 

.2000 

.  9838 

.9803 

.36 

.2500 

1.0272 

1.0037 

2.34 

.3000 

1.0491 

1.0425 

.64 

.3500 

1.0635 

1.0979 

-3.13 

.4000 

1.1015 

1.1736 

-6.15 

.4500 

1.1994 

1.2754 

-5.96 

.5000 

1.3866 

1.4112 

-1.74 

.5500 

1.6656 

1.5910 

4.69 

1 

! 


TABLE  VI 

Kj  (Nondimensional)  for  Three-Point  Bending,  S/W=8 


Crack  Length 

a/W 

Kj  (Nondimensional) 

Calculated  Using 
the  Weight  Function 

Kj.  (Nondimensional) 

Ref  28 

Relative 
Difference 
in  % 

.1000 

1.0630 

1.0156 

4.66 

.1500 

.9911 

1.0080 

-1.68 

.2000 

1.0314 

1.0183 

1.28 

.2500 

1.0755 

1.0438 

3.04 

.3000 

1.0964 

1.0839 

1.15 

.3500 

1.1087 

1.1403 

-2.77 

.4000 

1.1450 

1.2167 

-5.89 

.4500 

1.2431 

1.3189 

-5.75 

.5000 

1.4328 

1.4549 

-1.52 

5500 


1.7160 


1.6349 


4.96 


Fig.  5.11.  Stress  Intensity  Factor  (Nondimens ional)  for 
Three-Point  Bending,  S/W=8 
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VI.  Conclusions 


1.  This  work  shows  clearly  that  the  experimental 
weight  function  procedure  leads  to  excellent  results  for  the 
stress  intensity  factor  without  actually  testing  the  loading 
configuration  in  question. 

2.  A  set  of  experiments  of  a  single  loading  con¬ 
figuration  can  be  used  to  calculate  the  stress  intensity 
factor  solution  for  any  other  loading  configuration  applied 
to  the  same  geometry.  The  only  additional  information 
required  is  the  stress  distribution  at  the  crack  location 
for  a  noncracked  body. 

3.  One  should  limit  the  results  obtained  for  a 
different  loading  only  to  the  same  a/w  region  studied 

in  the  original  experiment.  Extrapolations ,  or  even  calcu¬ 
lations,  near  the  test  a/W  limits  may  lead  to  inaccurate 
results. 

4.  Because  data  scatter  can  lead  to  errors  when 
differentiated,  care  should  be  taken  to  collect  enough 
data  to  yield  reasonable  interpolating  functions. 

5.  Even  though  some  error  was  introduced  in  the 
Kj  and  tio  measurements,  the  numerical  procedure  using  the 
derivatives  of  Kj  and  n0  did  not  distort  the  Case  2  results 
beyond  the  basic  accuracy  of  the  Case  I  measurements. 


6.  The  laser  interferometric  technique  provides 
an  efficient  procedure  for  crack  opening  measurements  and 
leads  to  accurate  results  for  this  type  of  measurement. 

7.  The  Orange  (Ref  11)  conic  section  again  pro¬ 
vided  an  accurate  representation  for  the  crack  opening  pro¬ 
file  for  edge  cracked  strips.  This  conic  section  should, 
however,  be  compared  with  actual  profile  measurements  (as 
done  here)  prior  to  applying  it  to  other  geometries. 

8.  The  fact  that  the  actual  cracks  were  not  per¬ 
fectly  straight  on  a  micro-scale,  as  assumed  in  this  theory, 
is  the  most  likely  source  for  Kj  calibration  errors  for 
Case  1  loading  (four-point  bending) .  Perhaps  other 
materials  would  give  "straighter"  cracks  which  more  closely 
match  the  analytical  behavior. 
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Appendix  A 


The  Weight  Function 

The  idea  of  the  weight  function  has  been  intro¬ 
duced  by  Bueckner  (Ref  2)  and  was  also  discussed  in  (Refs 
2-8)  . 

The  following  analysis  shows  that  if  the  complete 
solution  of  the  stress  intensity  factor  and  the  crack  open¬ 
ing  displacement  for  a  crack  subjected  to  a  certain  load¬ 
ing  system  is  known,  then  the  stress  intensity  factor 
solution  for  the  geometry  under  another  loading  configura¬ 
tion  may  be  obtained  directly  from  the  known  solution. 

Let  us  consider  a  cracked  body  with  loads  P1,...,Pn 
as  described  in  Fig.  A.l. 


Fig.  A.l.  Loaded  Cracked  Body 
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The  elastic  energy  release  rate,  often  referred  to  as  the 
Griffith  energy  rate,  can  be  defined  following  Irwin 
(Ref  22)  as 


(A-l) 


where  U  is  the  total  potential  energy  stored  in  the  cracked 
body  and  a  is  the  crack  length. 

§  has  units  of  force  and  is  often  referred  to  as 
the  crack  extension  force 


U 


i  N 

i  I 


N 

I 


i=l  j=l 


C.  .P.P. 
l)  i  ] 


(A- 2) 


where  is  the  compliance  coefficient  defined  as  the 
deflection  at  point  i  due  to  a  unit  force  at  point  j . 


.  N  N  3C .  .  ,  N  9u, 

1  y  y  — il  P  P  =  I  y  — i 

2  i=i  j=i  9a  1  j  2i=i  9a 


p. 

X 

(A- 3) 


where 


N  N 

"  I  ui  =  l  CijPj 

1  j=l  1  j=l  -> 


(A-4) 


n?  is  the  displacement  at  i  due  to  load  Pj. 

Since  the  elastic  energy  release  rate  0  (based  on 
strain  energy  consideration)  and  the  stress  intensity  fac¬ 
tor  K  (based  on  crack  tip  stress  considerations)  are  both 
parameters  that  characterize  crack  growth,  a  relation  can 
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be  established  between  them  as  shown  by  Hertzberg  (Ref  15) . 


(A- 5) 


where  H  is  defined  as 


E  for  plane  stress 
E 

-  for  plane  strain 

1-v2 


(A-6) 


Since  K  is  the  linear  elastic  stress  intensity  fac¬ 
tor,  superposition  can  be  used  to  write 


N  N 

K  =  l  K.  =  l  k.  (a) P .  (A-7) 

i=l  1  i=l  1  1 


k^  stress  intensity  factor  per  unit  load.  Substituting, 
into  the  elastic  energy  release  rate 


.  N  N 

5  l  l 

H  i=l  j-1 


k±  (a)  kj  (ajP^ 


(A- 8) 


By  equating  the  Eq  (A-8)  to  Eq  (A-3)  we  conclude  that  they 
should  be  equal  term  by  term. 

Thus, 


kA  (a)  k,  (a) 
H 


1 

2 


5Cii(a) 

9a 


(A- 9) 


If  the  full  solution  for  K  is  known  only  for  one 
load,  say  P  ,  then 
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(A-10) 


k.(a) 


9C.  (a) 
1m 

3a 


1 

k  (a) 
in 


recall,  that  we  defined  in  Eq  (A-7)  k 

m 

and  in  Eq  (A-4)  C. 

un 


m 

1 

P 


m 


Therefore, 


and  since 


Then 


ki(a) 


H 

2 


~  m 
3ui 
3  a 


N 


K  =  l  k. (a) P . 
i=l  1  1 


K  = 


H 

2K 


n  3U; 


m 


y  — i —  p 
9a  i 


m  i=l 


(A-ll) 


(A-12) 


(A-13) 


For  an  arbitrary  distribution  of  P^,  and  a  set  of  surface 
tractions  T(s),  the  expression  for  K  can  be  written  as 


K 


T (s) ds 


If  one  defines  the  weight  function  as 


(A-14) 


hm(s,a) 


H  5um(s,a) 

2K  (a)  3a 
m 


(A-15) 


Then  we  can  readily  see  that  the  weight  function  depends 
only  on  the  geometry  of  the  cracked  body  (including  the 
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crack  length)  and  the  loading  configuration  denoted  by 
the  superscript  m. 

Substituting  the  term  for  the  weight  function 
into  Eq  (A-14) ,  we  get 


K 


hjn  ( s,a) T  (s)  ds 


(A- 16) 


We  can  therefore  conclude  that  the  stress  intensity  factor 
for  the  arbitrary  surface  traction  T(s)  can  be  obtained 
from  the  known  weight  function  we  computed  considering  the 
loading  configuration  m.  In  a  little  bit  more  general 
derivation,  the  body  forces  can  be  included  and  one  can 
write 


K 


T.hdr 


f  *h*dA 


(A-17) 


where  h  is  the  weight  function  obtained  from  any  known  load 
configuration , 

T  are  the  surface  traction  for  which  K  i;_  to  be 
determined, 

f  are  the  body  forces  for  which  K  is  to  be  determined, 

T  is  a  chosen  path  around  our  specimen  as  to  include 
all  the  surface  tractions  T,  and 

A  is  the  region  where  the  body  forces  act. 


61 


Appendix  B 

Justification  of  Assumptions 

Use  of  Linear  Elastic 
Fracture  Mechanics 

Linear  elastic  fracture  mechanics  (LEFM)  requires 
that  the  amount  of  plasticity  near  the  crack  tip  will  be 
relatively  small.  As  a  "rule  of  thumb"  LEFM  can  be  con¬ 
sidered  to  be  valid  if 

r  <  0.1a 
P  ~ 

where  r^  is  the  radius  of  the  plastic  zone  around  the  crack 
tip  and  a  is  the  crack  length.  Using  Irwin's  circular  model 
for  the  plastic  zone  incorporating  the  upper  bound  dimen¬ 
sion  as  obtained  by  a  plane  stress  solution  (see  Fig.  B.l), 


Fig.  B.l.  Irwin  Circular  Plastic  Zone  Model 
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r 


P 


_1_ 

2tt 


(B-l) 


•where  ay.s  t^ie  yi®ld  strength  of  the  material. 

Using  7075  T651 

ay*s  =  80.8  ksi 

The  maximum  used  in  our  tests  was  lower  than  15  ksi/Tn. 
Therefore,  the  maximum  radius  of  plastic  zone  that  has  been 
developed  was 


r 


P 


0.0055  in 


The  smallest  crack  length  used  was 


a  =  0.0815” 


Thus, 


'!*) 

^  a  / 


max 


0.0674  <  0.1 


The  conclusion  is  that  LEFM  was  valid  even  in  the  most 
extreme  conditions.  In  most  cases  the  crack  length  was 
significantly  longer  and  therefore  r  /a  smaller. 

r 


Plane  Stress/Strain  Considerations 

As  a  "rule  of  thumb"  (Refs  15-16) ,  if 


b  >  2.5 


(B-2) 
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where  b  is  the  material  thickness,  KIC  the  fracture  tough¬ 
ness,  and  a  the  yield  strength. 

y  •  S 

The  specimen  can  be  considered  thick  enough  to  be 
under  a  plane  strain  stress  field. 

for  7075  T651  plate  is  21. 4  ksi/Tn. 


Since  .24"-  .25" thick  specimens  were  used,  they  can 
be  viewed  to  be  under  plane  strain  conditions. 
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Appendix  C 

Laser  Interferometer  Data  Reduction 


A  typical  strip  chart  recorder  trace  of  fringe 
motion  crack  mouth  opening  and  load  is  shown  in  Pig.  C.l. 
From  this  chart  the  fringe  order  versus  load  data  can  be 
obtained  for  the  left  and  right  fringe  pattern  as  pre¬ 
sented  in  Table  C.l.  The  data  can  be  plotted  as  shown 
in  Pig.  C2. 

From  Fig.  C.2  slopes  of  Sm/P  of  the  left  and  right 
fringe  pattern  can  be  obtained  by  considering  only  the 
linear  part  and  ignoring  the  initial  nonlinear  part  of  the 
plot  that  is  associated  with  fatigue  crack  closure  effects. 
Those  slopes  are  averaged  to  eliminate  free  body  motion. 

The  experiment  is  repeated  and  results  are  averaged. 
The  value  of  the  slope  is  used  in  Eq  (4-1)  to  get  displace¬ 
ment  per  unit  load  data.  Through  Eq  (4-2)  stress  intensity 
factors  can  be  calculated. 


65 


A 


67 


Typical  Load-Fringe  Motion  Curve 
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Appendix  D 

Computer  Program  Description 

The  computer  program  has  been  written  to  perform 
the  numerical  computation  involved  in  the  process  of  data 
interpolation,  weight  function  calculations  and  finally 
evaluation  of  the  integral  in  Eq  (2-30)  to  determine 
for  Case  2. 

In  order  to  allow  maximum  flexibility  in  the 
usage  of  the  program  it  was  constructed  of  a  main  program 
that  performs  the  integration  and  nine  subprograms  that 
perform  all  the  other  necessary  calculations.  In  addition, 
two  library  subprograms  are  used,  one  from  the  IMSL  library 
and  the  other  from  the  CC  6600  library.  This  structure 
of  the  program  allows  using  it  for  different  cases  and  con¬ 
figurations  and  permits  the  user  to  enter  his  own  interpo¬ 
lation  schemes,  comparison  functions  and  the  noncracked 
body  stress  distribution  p(t)  that  makes  it  possible  to 
compute  any  desired  Case  2  stress  intensity  factor. 

The  program  is  written  in  FORTRAN  and  was  exe¬ 
cuted  on  the  CDC  system  at  Wright-Patterson  Air  Force 
Base.  Execution  time  for  the  three -point  bending  and 
uniform  tension  Case  2  stress  intensity  factor  varied 
from  two  to  four  seconds. 
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The  Main  Program 

The  main  program  reads  in  the  input  data,  calls 
in  the  different  subroutines  to  perform  necessary  computa¬ 
tion  and  performs  the  numeric  integration  by  using  the  IMSL 
library  subprogram  DCADRE.  An  iterative  scheme  "pushes" 
down  the  lower  integration  limit  until  the  desired  accuracy 
for  the  end  result  is  met.  At  the  end,  it  prints  out  com¬ 
puted  Kj  for  Case  2  compared  to  user-supplied  comparison 
data. 

Function  Subprogram  F (x) 

This  function  computes  the  value  of  the  inte¬ 
grant  (the  weight  function  multiplied  by  P(t))  at  each  dis¬ 
tance  t  from  the  crack  tip. 

Function  Subprogram  P (x) 

This  function  computes  the  value  of  the  stress 
at  the  crack  location  P(t)  for  a  noncracked  body  at  any 
distance  t  from  the  corresponding  crack  tip. 

The  function  should  be  supplied  by  the  user  and 
is  the  only  subprogram  that  has  to  be  altered  when  com¬ 
puting  Kj.  for  different  Case  2  configurations.  Care 
should  be  taken  as  related  to  the  form  in  which  P(t)  is 
supplied.  This  form  will  determine  the  form  of  the  calcu¬ 
lated  Kj.  A  factor  of  l/Zifa  is  already  included  in  the 
weight  function  calculation.  So  if  P(t)  is  supplied 
iirectly  in  terms  of  stress,  the  resulting  Kx  values  will 
■  *»  in  the  dimensional  form  of  Kj/Zra. 
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Library  Function  PLSCF 

A  CC  6600  Library  subroutine,  computes  least 
square  polynomial  fits  for  a  given  set  of  data.  Poly¬ 
nomials  up  to  6  degrees  or  Chebyshev  polynomials  of  any 
degree  can  be  specified  by  the  user.  For  further  informa¬ 
tion  see  (Ref  25) . 

Library  Function  Subprogram  DCADRE 

An  International  Mathematics  and  Statistics  IMSL 
Library  subprogram.  Performs  numeric  integration  of  a  func¬ 
tion  using  the  Romberg  extrapolation  scheme.  The  function 
DCADRE  is  capable  of  integrating  functions  with  "jump" 
discontinuities  and  certain  types  of  singularity.  This 
feature  is  most  important  in  our  case  because  of  the  crack 
tip  singularity  of  the  weight  function.  Even  if  the  desired 
accuracy  is  not  found,  DCADRE  returns  the  best  available 
estimate.  See  IMSL  Manual  (Ref  26)  for  further  information. 

Subroutine  KPOL 

This  subroutine  computes  a  fourth  order  least 
square  polynomial  fit  for  the  nondimensional  K^.  data  using 
PLSCF  subroutine  from  the  CC  6600  library.  It  prints  out 
the  polynomial  coef f icients C  ( 5). . . ,C (1)  and  returns  these 
values  to  the  main  program.  It  also  prints  out  the  error 
flag  IER  from  PLSCF. 

Subroutine  EPOL 

This  subroutine  computes  a  fourth  order  least 
square  polynomial  fit  for  the  nondimensional  nQ  data  using 
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PLSCF  subroutine  from  the  CC  6600  library.  It  prints  out 


the  polynomial  coefficients  D (5) , . . . ,D (1)  and  returns  them 
to  the  main  program.  It  also  prints  out  the  error  flag 
IER  from  PLSCF. 


Subroutine  EVEK 

Evaluates  the  nondimensional  Kj  from  the  fourth 
order  least  square  polynomial  and  the  derivative  of  Kj  with 
respect  to  crack  length  SKj/Sa  at  the  points  where  for 
Case  2  should  be  evaluated.  One  should  realize  that  the 
differentiation  considers  the  fact  that  the  nondimensional 
K-j.  usually  factors  out  l//a~  that  should  be  included  in  the 
differentiated  term.  The  derivative  3KI/3a  requires  the 
expression 


where  Y  is  the  nondimensional  KI  value  and  o  is  the  maxi¬ 
mum  tensile  stress  in  the  specimen.  and  SK^/Sa  are 

returned  to  the  main  program  in  a  dimensional  form  factored 
only  by  l/o  as  to  match  the  nQ  data. 


Subroutine  EVETTA 

Evaluates  the  nondimensional  n  from  the  least 
square  fourth  order  polynomial  and  the  derivative  of  no 
with  respect  to  crack  length  3f|Q/3a,  at  the  points  where 
Kj.  for  Case  2  should  be  evaluated.  One  should  realize 
that  the  differentiation  considers  the  fact  that  the 


72 


nondimens ional  n  is  factored  by  a  term  of  "a"  that  should 
o 

be  included  in  the  differentiated  term  as  applied  in  the 
above  paragraph. 

n  and  3n  /3a  are  returned  in  a  dimensional  form 
o  o 

factored  by  l/o  to  match  the  data. 

Subroutine  COPK 

Compares  the  nondimens ional  input  data  with  the 

interpolated  data  by  KPOL  at  the  test  points.  A  comparison 
table  that  includes  a  relative  error  in  percents  is  printed 
out  for  the  user's  evaluation  and  judgement. 

Subroutine  COPETA 

Compares  the  nondimens ional  nQ  input  data  with 
interpolated  data  by  EPOL  at  the  test  points.  A  comparison 
table  that  includes  a  relative  error  in  percents  is  printed 
out  for  the  user's  evaluation  and  judgement. 

Subroutine  KANA 

Evaluates  Kj  for  the  Case  2  loading  configuration 
from  a  user-supplied  comparison  function.  The  evaluation 
is  done  at  the  same  points  where  Case  2  is  evaluated 
from  the  experimental  data.  This  subroutine  can  be 
eliminated  if  the  above  comparison  data  is  available 
explicitly.  In  this  case  the  comparison  values  should  be 
assigned  directly  to  AK  array.  If  no  comparison  values 
are  available,  the  final  printout  table  should  be  altered 
to  avoid  undefined  terms. 
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List  of  Symbols  Used  in  the  Program 
(Alphabetical  Order) 

A  Constant  used  in  subroutine  KANA 

AERR  Input,  absolute  error  in  the  integral  routine 


AK (99) 


AKS  (99) 

B 

C  (5) 

D(5) 

DB 

DC 

DY (99) 
DYA(99) 

D2 (99) 
DZA (99) 

EDS 

ES (99) 


DCADRE 

K_  loA Ta  comparison  value  of  the  KT  Case  2 

r©  1 

nondimensional  evaluated  from  the  user-supplied 

comparison  function  at  points  defined  by  S 

Kt  /  a/i ra  calculated  values  of  Kx  Case  2  (non- 

*©  1 

dimensional)  at  points  defined  by  S 

Interim  term  for  weight  function  calculations 

Interpolating  polynomial  coefficients  for  KT 
data 

Interpolating  polynomial  coefficients  for  nQ  data 

Interim  term  for  weight  function  calculation 

3R/3a  interim  term  for  weight  function  calcu¬ 
lation 

Interim  derivative  term  for  DYA 

1/a  (SKj/Sa)  derivative  w-r-t  the  crack  length 
of  Kj.  Case  1  at  points  defined  by  S 

Interim  derivative  term  for  DZA 

1/a  (3n  /3a)  derivative  w*r-t  the  crack  length 
o 

of  riQ  at  points  defined  by  S 

Input,  desired  absolute  accuracy  in  KT  computa¬ 
tion 

(AKS-AK) 100/AK  relative  error  in  the  computed 


EST  Current  estimate  of  value  of  the  integral  by 

DCADRE 


EY(99)  (YC-Y) 100/Y  relative  error  in  percents  of 

interpolated  Kj  data  at  points  defined  by  S 
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EZ  (99) 


F 

G 

H 

IER 

M 


(ZC-Z) 100/Z  relative  error  in  percents  of 
interpolated  nQ  data  at  points  defined  by  S 

h,  the  weight  function  (subfunction) 

Upper  limit  of  the  integral 

Input,  H  material  property 

Error  flag  for  library  subprograms  PLSCF  and 
DCADRE 

Input,  number  of  experimental  data  points 


N 


Input,  number  of  points  where  K ^  should  be 
evaluated  for  Case  2 


P  Stress  distribution  for  the  noncracked  body  at 

the  crack  location  for  Case  2  loading  configura¬ 
tion  (subfunction) 

PI  ir 


Q 

R 

RELERR 


S  (99) 

T 

V 


Lower  limit  of  integral 

Interim  term  for  weight  function  calculation 

Input,  relative  error  in  the  integration  routine 
DCADRE 

Input,  a/w,  nondimens ional  crack  length  at 
point  where  experiments  were  performed 

Interim  term  for  weight  function  calculations 

Previous  interim  value  for  EST 


VA  (EST-V) ,  difference  between  last  and  current 

integration  estimate 

W(99)  Weight  corresponding  to  each  data  point  used 

in  PLSCF  (currently  defined  as  -1.  to  give 
equal  weight  to  each  point) 

WORK (25)  Work  array  for  PLSCF  subroutine 

X(99)  Input,  a/w,  nondimensional  crack  length  for 

points  where  Case  2  should  be  evaluated 

Y(99)  Input,  Kt/  a/ira,  nondimensional  KT,  experimental 

data  1 
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YA (99) 


KI/a,  interpolated  Kj  for  Case  1  at  points 
defined  by  x 

YC(99)  Kj/  a/ia  ,  interpolated  nondimensional  Kj  for 

Case  1  at  points  defined  by  S 

YK(99)  Kj/  a /Fa  ,  interpolated  nondimensional  Kj.  for 

Case  1  at  points  defined  by  x 

Z ( 99 )  Input,  n  *H/(2aa),  nondimensional  crack  mouth 

opening 

ZA ( 99)  no/a •  interpolated  nQ  at  the  points  defined  by  x 

ZC(99)  noH/(2aa),  interpolated  nondimensional  n  at 

points  defined  by  S 

ZK ( 99)  noH/(2aa),  interpolated  nondimensional  n  at 

points  defined  by  x 


Sequence  of  Input  Data 


Card 

Nunber 

Variable 

No.  of 
Variables 
per  card 

Real/ 

Integer 

Format 

1 

H 

1 

Real 

Unformatted 

2 

AERR,RELERR,EPS 

3 

Real 

Unformatted 

3 

M 

1 

Integer 

Unformatted 

3  to  (JVH-3) 

S (I) ,Y (I) ,Z(I) 

3 

Real 

Unformatted 

NH-4 

N 

1 

Integer 

Unformatted 

(M+5)  to  (MW+5) 

X(I) 

1 

Real 

Unformatted 

Sequence  of  Output  Data 

1.  Interpolation  polynomial  coefficients  for  K-j. 

2.  Interpolation  polynomial  coefficients  for  ri 

3.  Comparison  table  of  actual  to  interpolated 

4.  Comparison  table  of  actual  to  interpolated  nQ 
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5.  The  following  printout  will  appear  only  if 
singularities  have  been  met  by  subroutine  DCADRE.  Each 
time  a  singularity  was  met  a  warning  statement  is  printed 
out  followed  by  the  values  of  the  crack  length  (X(I))  the 
distance  from  crack  tip  (Q)  and  the  value  of  the  error  flag 
(IER)  for  which  the  singularity  occurred.  Values  of  XER 
equal  to  65  or  66  indicate  that  the  singularity  was  suc¬ 
cessfully  handled.  If  any  other  value  appears  consult 

the  IMSL  Manual  (Ref  26). 

6.  Comparison  table  of  for  Case  2  calculated 
values  and  supplied  comparison  values. 

Listing  and  Result  Sample 

The  following  pages  include  a  complete  listing  of 
the  program  and  a  sample  printout  of  results  for  uniform 
tension  as  Case  2  loading  configuration. 
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t  t 

5S 

H  W 
*-  X 

s* 

u  u 


a  vi 
x  UJ 

O  3 

o  -I 


O  3 
m  l* 

ar  z 

UJ  M 

§" 

z  ►. 


X  K 
o« 

«v  o 

U.  M 


n  r 
o  a 
«««  > 


or 

V)  o 
uj  ee 
m*  or 
h  jw 

M  * 

a  3  uj 

35£ 

3r* 

Z  JO 
m«»m 
VI  X  X 
IMS 


r  v>  *-  < 


•  «u 

•  • 
X  * 
•  -* 


•  XX 
X  •  - 

•  MV 


X  * 
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CO  HHO  N/OA  N/C • T  * 8,  R , 08 , OC 
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Crack  Shape — Conical  Approximation 


Following  Orange  (Ref  11) ,  the  shape  of  a  crack 
can  be  approximated  by  a  conic  section  equation.  Requiring 
that  the  crack  opening  at  the  crack  mouth  will  match  the 
actual  crack  mouth  opening  no  one  can  write  a  general  form 
for  the  conic  section;  equation; 


(E-l) 


✓ 


Fig.  E.l.  Crack  Opening  Conical  Approximation 

where  t  is  the  distance  from  the  crack  tip,  a  is  the  crack 
length,  and  m  is  the  conic  section  coefficient  to  be  deter¬ 
mined.  The  type  of  conical  section  depends  on  m  if 


m  <  0 


ellipse 
m  =  0  -  parabola 

0<m<°»  -  hyperbola 

m  =  co  -  a  pair  of  straight  lines 

In  order  to  determine  the  unknown  conic  section 
coefficient  in,  let's  match  the  radius  of  curvature  at  the 
crack  tip  from  the  Eq  (E-l)  to  the  curvature  resulting 
from  the  crack  tip  displacement  field  as  expressed  in  Eq 
(2-16)  . 

Using  the  following  definition  for  the  curvature 
of  a  two-dimensional  curve  defined  by 


y  =  f(x) 

the  curvature  K  is  given  by 


d2x 

dy2 


and  the  radius  of  curvature  is 


(E-2) 


(E-3) 


K  -  [|] 
provided  K  ^  0 


If  one  uses  the  near  crack  tip  displacement  field  from 


t 


3t  _  ttH2 

an  4Kj 


n 


and  at  the  crack  tip 


n  =  o 
t  =  o. 


(E-5) 


(E-6) 


Prom  Eq  E-6) 


3t 

3n 


n=o 


o 


(E-7) 


Also  from  Eq  (E-6) 

lit  =  mi  (E_8) 

3n2  4Ki 

Substituting  Eqs  (E-7,  E-8)  into  Eq  (E-3)  and  then  using 
Eq  (E-4) ,  the  radius  of  curvature  at  the  crack  tip  from 
the  crack  tip  displacement  field  is 


R 


(E-9) 


Now  find  the  radius  of  curvature  from  the  conic  section 
Eq  (E-l) 

2n  2  mn^2 

n2  =  - ° —  t  +  — ° -  t2  (E-10) 

(2+m)  a  (2+ifi)  a2 
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Differentiate  twice  with  respect  to  n 


2n 


2n  = 


(2+m)  a 


a.  2itin  2 
9t  ,  O 

9h 


( 2+in)  a2 


.  9 1 

t  *  9^ 


(E-ll) 


2n. 


9  2 1 


2  = 


(2+m) a  3n2 


(2+m)  a2 


[(#  •  £?] 


(E-12) 


From  Eqs  (E-ll,  E-12)  at  the  crack  tip 


t=0,  n=0 


The  value 


and 


9_t 

9n 


0 


3 2 1  _  (2+m) a 
3n2  ti  2 


(E-13) 


If  Eqs  (E-3,  E-4)  are  used,  the  radius  of  curvature  at  the 
tip  of  the  conic  section  is 

n  2 

R  =  —  (E-14) 

(m+2) a 


Now,  matching  the  radius  of  curvature  at  the  crack  tip  of 
the  conic  section  (Eq  E-14)  to  the  one  that  is  derived  from 
the  crack  tip  displacement  field  Eq  (E-9) ,  one  can  determine 
the  conic  section  coefficient 


m  = 


2 


(E-14) 
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Thus,  following  Orange  and  substituting  the  stress  inten¬ 
sity  factor  Kj,  with  the  nondimensional  stress  intensity 
factor  Y  defined  by  the  relation 


KI  =  Ya/a 
The  value  of  in  becomes; 


2 


(t  15) 


where 


Y  = 


K 

a/a 


(E-16) 
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